Regular Spectra and Universal Directionality of Emitted Radiation from a 

Quadrupolar Deformed Microcavity 



Jeong-Bo Shim and Hai-Woong Leqj 

Department of Physics, Korea Advanced Institute of Science and Technology, Daejeon, 305-701, Korea 



o 
o 

(N 



(N 



C/3 

o 

> 

q 

O 

Oh' 



m 
o 

o 

o , 

• I— I 

CZ3 ■ 



X 



Sang-Bum Lee, Juhee Yang, Songki Moon, Jai-Hyung Lee, and Kyungwon An 
School of Physics, Seoul National University, Seoul, 151-742, Korea 

Sang Wook Kim 

Department of Physics Education, Pusan National University, Busan 609-735, Korea 

(Dated: February 2, 2008) 

We have investigated quasi-eigenmodes of a quadrupolar deformed microcavity by extensive nu- 
merical calculations. The spectral structure is found to be quite regular, which can be explained 
on the basis of the fact that the microcavity is an open system. The far-field emission directions 
of the modes show unexpected similarity irrespective of their distinct shapes in phase space. This 
universal directionality is ascribed to the influence from the geometry of the unstable manifolds in 
the corresponding ray dynamics. 



I. INTRODUCTION 

Symmetrical (cylindrical or spherical) microcavities 
have attracted much attention in the past as laser res- 
onators, due to their compact and simple geometry, eas- 
iness to fabricate, and ultra-high Q values, which is at- 
tributed to formation of the so-called whispering gallery 
modes (WGM's) 0. The isotropic radiation pattern of 
light emitted from the symmetrical microcavity, however, 
reduces their practical usefulness. One way to overcome 
this problem, as has been proposed by Nockel et al., is to 
deform the shape to construct asymmetric resonant cavi- 
ties (Arc's) i2fci2ij- Directional emission has indeed been 
observed experimentally from asymmetrically deformed 
microlasers made of semiconductors 0, |^, dye jets @ 
and polymers 7]. 

Fundamental optical properties of ARC's can be partly 
understood in terms of the ray dynamics. From the clas- 
sical dynamical viewpoint, it belongs to a class of sys- 
tems having mixed-type phase space, i.e., a mixture of 
both regular and irregular trajectories H,0j0|- When 
deformation is very small, the phase space is dominated 
by regular rays of whispering gallery type and these rays 
escape from the cavity only by evanescent leakage (tun- 
nelling). At larger deformations, however, emission prop- 
erties are mainly determined by chaotic rays that diffuse 
stochastically and refract out of the resonator when the 
angle of incidence x reaches the critical angle Xc d-e-j 
sinxc — ^/n, where n is the refractive index of the res- 
onator. This simple ray dynamics model predicts that at 
large deformations the escape occurs primarily into the 
tangential direction near the points of maximum curva- 
ture 2, 11]. 

Although the model captures the essence of the direc- 
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tional property of emitted radiation from ARC's, there 
may exist other complications. For example, a phe- 
nomenon referred to as dynamical eclipsing |l2l | occurs 
when stable islands occupy the phase space region at 
which the escape would have taken place without them, 
resulting in strong suppression of emission intensity in 
directions predicted by the ray model and appearance of 
split peaks in nearby directions. There are evidences that 
stable or unstable periodic orbits strongly influence the 
emission pattern of ARC. The emission originating from 
bow-tie modes 0, 0| , whose corresponding classical or- 
bit is a stable island, and the so-called scarred modes of 
hexagonal and triangular ||| unstable periodic orbits, 
has been observed. So far the direction of the far-field 
emission obtained in experiments has been a commonly 
used tool to investigate the characteristics of the lasing 
mode of ARC, since it has been assumed that the shape 
of the wavefunction inside the cavity strongly influences 
or sometimes even completely governs the far-field emis- 
sion pattern. It has been reported, however, that, due to 
the correction with Fresnel coefficients, the emission di- 
rection of the mode with a rather smaller Q value can be 
considerably deviated from an expectation based upon 
the above simple assumption Deviations from the 
prediction of the simple ray dynamics model also occur 
when the nonlinear effect of a medium is involved |l4j. 

The simple ray dynamics model has recently been re- 
fined by Schwefel et al. . They found that highly direc- 
tional emission patterns from ARC's persist well beyond 
the deformation limit predicted from the ray model 0, 
and developed a more accurate model which emphasizes 
the importance of unstable manifolds of short periodic 
orbits. However, they did not show the direct evidence 
of the crucial role played by unstable manifolds in con- 
nction with their structure embedded in quantum wave- 
functions. A similar analysis has been reported by Lee 
et al. in a stadium shaped microcavity, although the au- 
thors did not precisely mention the role of the unstable 



2 



manifolds [TH . 

The main purpose of this paper is to investigate opti- 
cal properties of a quadrupolar micro-resonator with low 
refractive index (n — 1.361) by extensive numerical cal- 
culations. The reason why we choose the specific refrac- 
tive index is that this work is initiated from our previous 
experiment of the dye-jet microcavity laser We would 
like to deliver two main results: firstly the energy level 
distribution shows deviations from an expectation from 
the usual random matrix theory (RMT), and secondly 
the far-filed emission direction exhibits a nontrivial uni- 
versal feature. 

To a certain degree, optical properties of ARC's can be 
understood by using an analogy with chaotic dynamical 
properties of the billiard system which has been exten- 
sively studied in the context of nonlinear dynamics. The 
quantum mechanical manifestation of classical chaos, re- 
ferred to as quantum chaos, has generally revealed that 
there is a close relation between the classical nonlinear 
dynamics (whether it is regular or chaotic) and quan- 
tum spectral statistics (whether the energy level spa cing 
distribution is described by Poissonian or Wigner) [l6j |. 
Based on this observation, one might expect that there 
also exists close correlation between the ray dynamics 
and the statistical properties of the eigemodes of ARC's. 
Specifically, energy levels may be regularly distributed at 
small deformations, where the ray dynamics is predom- 
inantly regular, while they may look quite irregular at 
large deformations where chaotic rays dominate. We will 
show, however, that this simple expectation does not ex- 
actly hold true. There is an essential difference between 
ARC and the billiard system: The former is an open 
system while the latter is a closed one. The openness 
of ARC raises a nontrivial question: Where is the emit- 
ted output directed? This is completely meaningless in 
a closed system. Although the analogy between ARC 
and the billiard system is a useful starting point because 
the billiard has been thoroughly studied in the past, one 
cannot push the analogy too far. 

This paper is organized in the following way: In 
Section II, we briefly explain both the ray dynamics 
and the wave properties of a quadrupolar billiard, and 
a quadrupolar deformed microcavity (QDM). In Sec- 
tion III, we show numerical results for spectra and far- 
field emission directionality of a QDM and discuss the 
regularity of the spectra and the universal directionality 
we obtained from the numerical computations. Finally 
in Section IV, we present a conclusion. 



II. A QUADRUPOLAR DEFORMED 
MICROCAVITY 

When a liquid jet with circular cross section is de- 
formed, its cross section can in general be described by 




FIG. 1: A schematic view of BirkhofT coordinates: They con- 
sist of the distance along the boundary S{(j)) and the reflection 
angle, x- To make them canonically conjugate, sinx is used 
instead of x- 

the following multipole expansion 

-j^ oo oo 

= + X! cos -I- ^ bj sin j0, (1) 

where r is the distance from the center and (p the 
azimuthal angle. The lowest order even-symmetric 
quadrupolar contribution becomes dominant when the 
perturbation inducing deformation is small. The cross 
section of the dye jet can then be described by 

r(^) = -=lL==(l+ecos2</.), (2) 

where vq and e are respectively the radius of the unde- 
formed (circular) cavity and the deformation parameter. 
The denominator in Eq. ||2J) guarantees that the cross 
sectional area is conserved irrespective of e, which is con- 
sistant with the experimental condition of Ref. @. In 
this paper we focus on such a specific ARC, i.e., a QDM. 

A. A quadrupolar billiard: closed system 

At first, we consider the ray dynamics of a QDM with 
hard walls, namely a quadurpolar billiard. The Birkhoff 
coordinate system is quite useful to study the ray dy- 
namics in a billiard. It consists of two variables S (or 
equivalently (/)) and sinx, which respectively represent 
the arc length of the position at which the ray strikes 
the boundary and the sine of the incident (or reflected) 
angle at the boundary (Fig. QJ. In Fig. |21we show 
the Poincare surfaces of section (PSOS) of the dynamics 
in the quadrupolar billiard by using Birkhoff coordinates. 
It exhibits a typical Kolmogorov-Arnold-Moser scenario 
[T^ . i.e. the evolution toward stochastic behavior as the 
deformation parameter e is increased. 
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FIG. 2: The Poincares surfaces of section of a quadrupolar billiard for (a) e = 0, (b) e — 0.05, (c) e = 0.10, and (d) e = 0.12. 



The wave nature of a quadrupolar billiard can be de- 
scribed by the following Hclmholtz equation with Dirich- 
let boundary condition 



V^^{r) + k^^Pir) = 0, (3) 



where k is the wave number. For a circular billiard, the 
solution of Eq. Intakes the form Jm(r) exp(«TO(/)), where 
Jm is a Bessel function. All the eigenvalues ki^ are ana- 
lytically determined and form a well organized structure 
since they are just the Ith zeros of the mth order Bessel 
functions. When a circular billiard is deformed, how- 
ever, the distribution of eigenvalues starts to deviate from 
the regular structure. For large enough deformations, in 
which case the corresponding ray dynamics is predom- 
inantly chaotic, the well established RMT plays an im- 
portant role. Although the nearest neighbor level spacing 
distribution of eigenvalues is described by the universal 
Wigner function, the eigenvalues themselves show quite a 
complex structure due to a lot of avoided crossings. Fig.|31 
shows the evolution of eigenvalues (fcr) as a function of 
the deformation parameter. 
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FIG. 3: Evolution of the eigenvalues as a function of e. It 
shows a complicated structure. 



B. A quadrupolar deformed microcavity: open 
system 

A QDM is an open system in the sense that the ray 
with the incident angle larger than a critical angle is al- 
lowed to escape from the cavity. It does not make any 
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FIG. 4: Evolution of the real part of the eigenvalues of the 
mode with I = 1 (the dotted line), I = 2 (the solid line), and 
/ = 3 (the dashed line) of a QDM as a function of e. 



influence on the ray dynamics itself inside the cavity, but 
reduces available phase space: e.g. for n = 1.361 the 
phase space area below sin^c = 1/n « 0.73 can be ig- 
nored. 

The wave properties of an open system somewhat dif- 
fer from those of a closed one. The eigenvalues of Eq. Q 
now take complex values. Their imaginary part is asso- 
ciated with the decay rate of the mode. To quantify the 
decay rate, a cavity quality factor Q defined as —2xkr/ki, 
where kr and ki are respectively real and imaginary parts 
of the eigenvalue, is commonly used . The real part of 
the eigenvalue is represented by the size parameter de- 
fined as nkrVQ, which is abbreviated as nkr if there is no 
confusion. 

As far as a circular microcavity is concerned, Eq. 
is still separable so that two good quantum numbers, ra- 
dial (I) and angular ones (m), can be well defined. The 
eigenvalues kim are analytically given by matching be- 
tween Bessel and Neumann functions. Similarly to a cir- 
cular billiard, one can see the well-defined free spectral 
range (FSR) among the modes with the identical radial 
quantum number /. 

Figure^presents a continuous variation of the real part 
of eigenvalues as the deformation is increased. Surpris- 
ingly, the eigenvalues vary quite regularly in comparison 
with the complicated structures in Fig. O This regu- 
lar structure remains, although the deformation becomes 
large enough to generate global chaos in classical dynam- 
ics. The monotonic decrease of nkr implies the red shift 
of spectral lines. In case of WGM it can be explained 
by considering an elongated optical path due to defor- 
mation. In order to explain the regularity of spectrum, 
we need to investigate more carefully the nature of the 
modes. In the following section their phase space distri- 
bution will be discussed based mainly upon the Husimi 
plot. 



When the system is open, a question naturally arises: 
Into which direction is the output field emitted? In a 
circular cavity, the emission direction should be isotropic 
due to the circular symmetry, while it is highly non-trivial 
to address it for ARC's. In the viewpoint of the ray 
dynamics, the emission direction is determined from the 
position (0,sinx) at which the incident angle of the ray 
becomes equal to the critical angle. In wave mechanics, 
however, one should find quasi-eigenmodes, which are an 
open system analogue to eigenmodes of a closed system, 
by direct numerical calculation. In this work we use the 
boundary element method to compute quasi-eigenmodes 

MM. 



III. REGULAR SPECTRUM AND UNIVERSAL 
EMISSION DIRECTIONALITY 

In this section we discuss the origin of the regularity of 
the spectrum and the universal emission directionality. 
The Husimi distribution function is a very powerful tool 
to study these issues. 



A. The Husimi plot of the eigenmode of the 
deformed microcavity 

The Husimi distribution function can be regarded as a 
mathematical description of the phase-space distribution 
in which the minimum-uncertainty Gaussian wave packet 
is used as a basis l23: 



H^iqa,ko) = 
where 

£.{q - 90, fco) = exp 



■0(9) C{q-qo,ko) dq 



(q - qo? 

2a 



ikoiq - qo) 



(4) 



(5) 



In order to obtain the Husimi plot in Birkhoff coordi- 
nates, it is necessary to modify Eqs. Q and (0). Taking 
it into account that is 27r periodic, the basis function 
^ can be rewritten as 



<^((/) - (j)Q,smx) = ^exp 



2a 

+iksmx{(t) - 0o)] 



(6) 



and the integration in Eq. should be taken from to 
27r [23, 13 • F'or our calculation wc set the width of the 
gaussian as cr = ^/2/ko. 

In Fig. |31 we present Husimi plots of several quasi- 
eigenmodes obtained numerically. The I = 1 modes in 
a circular cavity evolves from a horizontally flat shape to 
the distribution localized somewhere in the phase space 
in Fig. ISJc). A similar evolution is also observed in the 
I = 2 modes except that it is localized on the period-6 
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FIG. 5: Husimi plots of a QDM with (a)e = 0.05, (b) 0.1 and (c) 0.16 for the / = 1 modes, and with (d) e = 0.05, (d) 0.1 
and (f) 0.16 for the I = 2 modes. The corresponding complex size parameters are (a) 96.64 — iO.003, (b) 96.89 — z0.004, (c) 
95.87 -i0.004, (d) 96.10 -i0.005, (e) 96.25 - iO. 014, and (f) 96.61 - i0.029. The critical line defined by sinxc = 1/n is indicated 
with the red horizontal lines. 



orbit which loses its stabihty around £ ~ 0.12. The I = 2 
modes then forms a hexagonal shaped scarred mode ex- 
tensively investigated in our previous work 0. All the 
modes identified with Z = 2 in Fig. 0] appear to have 
the similar shapes in Husimi plots as shown in Fig. 
It is more difficult to identify the corresponding classical 
orbit of the I = 1 modes. It seems to correspond to the 
period-6 orbit based upon Fig. EJc), since it has six en- 
hanced localized probability maxima in the phase space. 



However, their locations along the vertical axis (sinx) fit 
themselves to those of the period-8 orbit. Interestingly, 
this mode exactly corresponds to the period-8 orbit for 
larger nkr. We believe that the I = 1 mode is experi- 
encing the transition from the period-6 to the period-8 
orbit. More details will be discussed elsewhere. 
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FIG. 6: Husimi plots of the I = 2 mode at (a) nkr — 96.073 and (b) nkr — 98.12 for e = 0.16. These two modes are separated 
by two units of FSR. 



B. Regular spectra 



C. Universal emission directionality 



In closed systems, chaotic states which are irregularly 
distributed over available phase space are frequently ob- 
served, so that they strongly interact with and repel each 
other to generate avoided crossings and form complicated 
spectral structures as shown in Fig.O In contrast, quasi- 
eigenmodes in open systems hardly form chaotic states. 
The reason is that the chaotic modes have too small Q 
values to form modes in the practical sense, since they de- 
cay in a short time due to refractive escape. In Fig.2]the 
I = 3 mode seems to disappear at e ~ 0.08, which implies 
that it is difficult to find the mode due to a very small Q, 
since, for a given nkr, the larger / is the smaller Q is. The 
mode with even larger I (> 3) may exist in Fig. O but is 
practically meaningless, because Q is extremely low. One 
may more frequently find chaotic modes for larger nkr, 
because Q is exponentially enhanced as nkr is increased. 
It leads us to the conclusion that in open systems chaotic 
quasi-eigenmodes are rarely observed unless nkr is large 
enough. 



It is no wonder that non-chaotic modes show devia- 
tions from the prediction based upon RMT, since they 
are localized somewhere in phase space. It is unlikely 
that for a given nkr the I = 1 mode has a similar shape 
with the / = 2 mode. In the semiclassical theory the 
larger / the bigger the incident angle of the ray. It implies 
that the modes with different Ps form distinct shapes in 
phase space. It is clearly shown in Figs.|3c) and (f) that 
those two modes hardly overlap with each other. As e 
is increased, they independently evolve with their own 
shape maintained in phase space. That explains the reg- 
ular spectral structure observed in Fig. ^ 



The quasi-eigenmodes of an open system contain 
evanescent waves outside the cavity, which cause decay 
of the modes in time. Figures C^a) and (b) show that 
some of the waves leak out, which cannot be seen in the 
Husimi plot. The output directionality is in practice the 
most important characteristics that we seek from ARC's. 
It is more definitely illustrated by the far-field emission 
intensity as a function of the output direction as shown in 
Figs. c) and (d) . Surprisingly, these two figures demon- 
strate an almost identical pattern whose maximum is lo- 
cated around 40°. It should be emphasized that the cor- 
responding Husimi plots of these two modes differ from 
each other as shown in Figs. \^c) and (f). Since it has 
been believed that the emission direction is determined 
from the internal shape of the quasi-eigenmode, it is nec- 
essary to explain this unexpected observation. 

The clue for the origin of the identical far-field direc- 
tionality irrespective of the shapes of a quasi-eigenmode 
comes out of the fact that the far-field pattern obtained 
from ray dynamics shown in Fig. IHl^d) is not much dif- 
ferent from those from the wave calculations. In the ray 
calculations, the ray incident with the angle smaller than 
the critical angle is reflected in or emitted out with the 
probability determined by Fresnel coefficients. 

In a long time the ray dynamics in the chaotic region is 
governed by the so-called unstable manifolds T^. They 
are easily seen by launching a set of localized initial con- 
ditions of rays in quadrupolar billiard and following the 
dynamics for a few tens of collisions [see Figs. Ela) and 
(b)]. As time goes to the infinity, the ray dynamics in a 
closed system will finally become ergodic so that the rays 
are uniformly distributed over phase space. As far as an 
open system is concerned, however, rays escape from a 
cavity before reaching uniform distribution. Therefore, 
the available phase space is limited, and the dynamics in 
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FIG. 7: Probability distributions of the wavefunctions of the 
quasi-eigenmodes in real space for (a) nkr = 97.72 {I = 1) and 
(b) nkr = 98.12 (Z = 2). The far-field emission intensities are 
shown for (c) the mode corresponding to (a) and (d) the mode 
corresponding to (b). 

the restricted space is described by the niotion along a 
few dominant unstable manifolds. As a result, the unsta- 
ble manifolds dominate the long (but not extremely long) 
time dynamics in ARC's, and consequently the output 
directionality. In Fig. IHl^c), we plot the accumulated in- 
tensity in phase space, obtained by recording the Birkhoff 
coordinates at which each ray escapes and totaling the 
number of the escaped rays weighted by Fresnel refrac- 
tive coefficients at each phase space point (j^^sinx). The 
intensitty pattern resembles the structure of the unstable 
manifolds shown in Fig. 8(b). In Fig. 8(d) we show the 
far field emission intensity obtained from the ray dynam- 
ics calculations, which is seen to be consistent with the 
results of the wave calculations shown in Figs. 7(c) and 
7(d). 

The unstable manifolds can also be found through a 
careful examination of Husimi plots, as Figs. El^a) and 
inib) indicate. Such a fine structure could not be seen 
in the original Husimi plots of Figs. El because it is as- 
sociated with a tiny probability. No matter how small 
the probability associated with the unstable manifold is, 
it significantly overlaps the region below the critical an- 
gle. It means that this fine structure mainly contributes 
to the output directionality regardless of the tiny prob- 
ability, unless the Q value is not too small. Note that 
Q ~ 10'' for both the / 1 and Z = 2 modes. We 
may now conclude that, based on both the ray dynamics 
calculations and the wave calculations, the emission di- 



rectionality is determined mainly by the structure of the 
unstable manifolds. 

In order to show that the above conclusion is quite 
general, we also investigate the Husimi plots and the far- 
field emission directionality of the modes for larger nfcr, 
namely ~ 200. FigureslTUlfaVfe) show a variety of shapes 
of Husimi plots of the modes: For example, (b) is an 
octagonal and (d) is a hexagonal scarred mode. Another 
interesting point found in Fig. (10) is that each mode 
does not have any considerable overlap with the others. 
In this sense, every mode in Fig. 1101 is quite distinct from 
one another. Surprisingly again, the far-field emission 
directions are almost identical for all such five distinct 
modes as shown in Fig. 1111 As mentioned above, it is 
also attributed to the geometry of the unstable manifolds 
as clearly seen in Figs. llon Vfi). Therefore, we call the 
identical far-field emission directionality determined from 
the unstable manifolds a universal directionality. 



IV. CONCLUSION 

In this work we have studied spectral properties and 
mode structure of a quadrupolarly deformed microcavity 
and the output directionality of the radiation emitted 
from it. We summerize the main conclusion of our work. 

(1) When the system is open like a QDM, chaotic or- 
bits no longer easily form a quasieigenmode since their 
Q values are expected to be very low due to the refrac- 
tive escape. Consequently the localized modes such as 
scarred ones are more frequently observed in comparison 
with the usual closed chaotic system. This explains the 
regular structure of the spectrum with well-defined level 
spacings even at large deformations. 

(2) The direction of the radiation emitted from a 
quadrupolar microcavity resonator is determined mainly 
by the geometry of the unstable manifolds of the res- 
onator. As the geometry of the unstable manifolds is 
governed by classical dynamical properties of the res- 
onator and thus is independent of the structure of the 
excited modes, the output directionality exhibits univer- 
sality, largely independent of the degree of deformation 
and the pump power. 
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refractive coefficients, (d) The far-field emission intensity obtained from the ray dynamics. 
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unstable manifolds are shown in the background with black curves. 
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